F/G  20/10 


CO 

CO 

a 

co 


Office  of  Naval  Research 


Contract  N0ppi4  78-C-.0724 


4  Ju 1*3*981 

Principal  Investigator 
J.  C.!  Rasaiah 

-h 

Insti tution 
University  of  Maine 
Orono,  Maine  04469 


s 


DTIC 

cC'iT 

SEP  9 1981 

H 


CJ> 

ixJ 

U.J 

d 


v  i 


ter  public  release; 
tion  Unlimited 


81  8  20  191 

-  /  7  ' 


j 


b.  Publication  Citations 


1.  The  adsorption  of  dipoles  at  a  wall  in  the  presence  of  an  electric 
field:  The  RLHNC  Approximation  Journal  of  Chemical  Physics  73, 
3980,  (1980).  John  M.  Eggebrecht,  Dennis  J.  Isbister  and 
Jayendran  C.  Rasaiah. 

2.  Electrostriction  and  Dipolar  Ordering  at  an  Electrified  Wall, 
Chemical  Physics  Letters  79,  189,  (1981).  J.  C.  Rasaiah, 

D.  J.  Isbister  and  G.  S tel  1 . 

3.  Nonlinear  Effects  in  Polar  Fluids.  A  Molecular  Theory  of  Electro¬ 
striction.  Journal  of  Chemical  Physics  (1981)  -  to  appear. 

J.  C.  Rasaiah,  D.  J.  Isbister  and  G.  Stell. 


c.  Data  on  Scientific  Collaborators 

1.  Mark  Cilley  -  Programmer 

2.  John  M.  Eggebrecht  -  Graduate  Student 

3.  Dennis  J.  Isbister  -  Postdoctoral  Research  Associate 

(from  Chemistry  Department,  Faculty  of  Military  Studies, 
University  of  New  South  Wales,  Duntroon,  ACT,  Australia) 

4.  Noel  E.  Thompson  -  Visiting  Professor 

(from  Department  of  Mechanical  Engineering,  University 
of  New  South  Wales,  Duntroon,  ACT,  Australia) 


w 


Summary 

Work  completed  under  this  grant  on  polar  fluids  includes  an  extension 
of  perturbation  theory  to  polar  molecules  with  point  octopoles  and 
hexadecapol es  as  well  as  a  preliminary  investigation  of  molecules  with 
non-spherical  cores.  In  addition,  the  ordering  of  dipolar  molecules  near 
a  flat  electrified  wall  has  been  investigated  in  the  linearized  hypernetted 
chain  (LHNC)  and  quadratic  hypernetted  chain  (QHNC)  approximations.  This 
has  led  to  a  new  study  of  el ectrostriction  and  non-linear  effects  in  the 
polarization  density  of  dipolar  fluids. 


This  report  describes  research  on  polar  fluids  carried  out  at  the 
University  of  Maine  at  Orono  under  contract  from  the  Office  of  Naval  Research 
Apart  from  the  principal  investigator  (Professor  J.  C.  Rasaiah),  a  visiting 
fellow.  Professor  Noel  Thompson  from  Australia,  a  graduate  student,  Mr.  John 
Eggebrecht  and  a  postdoctoral  fellow.  Dr.  Dennis  Isbister  from  the  Royal 
Military  College,  University  of  New  South  Wales  in  Australia,  worked  on  this 
project;  the  last  two  receiving  partial  financial  support  under  the  present 
contract. 

The  primary  objective  of  this  investigation  was  to  extend  the  theory 
of  polar  fluids  developed  by  Stell,  Rasaiah  and  Narang^  to  (a)  spherical 
molecules  with  point-octopoles  and  poi nt-hexadecapol es  which  could  serve 
as  a  reasonable  model  for  methane  and  sulfur  hexafluoride;  (b)  non-spherical 
molecules  with  dipoles  and/or  quadruples  such  as  HCN  and  0.,.  A  second 
objective  was  to  extend  our  study  of  dipoles  to  surface  adsorption  near 
an  electrified  wall  . 

A  Monte  Carlo  calculation  of  the  free  energy  changes  that  accompany 
the  switching  on  of  polar  interactions  between  the  molecules  was  planned 
as  an  unambiguous  test  of  the  accuracy  of  the  theory  when  applied  to  non- 
spherical  cores.  Some  aspects  of  this  work  are  still  in  progress,  while 
parts  (a)  and  (b)  of  the  research  described  have  been  completed.  In 
addition,  a  very  promising  and  important  theoretical  investigation  of  the 
adsorption  of  dipoles  at  a  wall  in  the  presence  of  an  electric  field  has 
been  initiated  with  the  assistance  of  Dr.  Dennis  Isbister  and  John 
Eggebrecht.  In  what  follows,  the  research  that  has  been  completed  with 
support  from  ONR  is  described  in  greater  detail  and  future  plans  to  continue 
this  work  are  outlined. 


-  2  - 


I.  Perturbation  Theory  for  Polar  Fluids 

While  the  potential  energies  of  interaction  amongst  point-dipoles 

and  point-quadrupolesare  well  known,  there  are  inconsistencies  and  ambigu- 

2 

ities  in  the  literature  when  higher  order  multipoles  were  included. 

It  was  therefore  decided  to  calculate  these  energies  from  first  principles, 
limiting  our  attention  only  to  octopolar  and  hexadecapolar  interactions 
beyond  quadrupolar  charge  distributions.  It  became  necessary  therefore  to 
perform  the  tedious  contractions  of  multipolar  and  field  tensors  the 
eight  index  field  tensor  for  the  hexadecapole-hexadecapole  energy  contains 
nearly  700  terms.2  The  two-body  potentials  beyond  the  simple  dipole-dipole, 
dipole-quadrupole  and  quadrupol e-quadrupole  interactions  were  found  to  be 
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where  p,  e,  a,  and  i  are  the  dipole,  quadrupole,  octopole  and  hexadecapole 
moments  and 
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where  (e..,(*>.)  are  the  polar  and  azimuthal  angles  of  orientation  of  the 
axis  of  symmetry  of  the  ith  molecule. 

The  potential  energy  functions  were  then  expanded  in  spherical 
harmonics,  to  allow  analytic  evaluation  of  the  angular  integrals  in  the 
x-expansion  of  the  free  energy 
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where  u.  =  ( 6 ^ ^ ,  and  S^1  (u^)  are  the  spherical  harmonics,  the 
spatial  and  angular  variables  of  the  potential  energy  can  be  separated. 

Using  the  orthogonality  properties  of  the  spherical  harmonics,  and  the 
previously  derived  potential  energy  functions  we  have  found  the  following 
expansion  coefficients  for  the  higher  order  multipole  interactions  (t-j  =  1, 
2,  3,  and  4  correspond  to  dipole,  quadrupole,  octopole  and  hexadecapole 
moments ) 
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The  general  pattern  of  these  coefficients  allows  a  check  on  the  tedious 
and  mistake-prone  calculations. 
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With  these  results  the  terms  of  0(x  )  in  the  free  energy  expansion  can 
be  calculated.  To  go  beyond  this,  using  a  Pade  approximant  for  the  free 
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energy  in  the  manner  proposed  by  Stell,  Rasaiah  and  Narang,  requires  at 
least  the  term  of  0(x  ).  To  this  end  the  corresponding  angularly  averaged 
three-body  potentials  for  higher  multipoles  were  needed.  They  were  not 

3 

available  in  the  literature  but  were  derived  in  the  manner  of  R.  J.  Bell 
and  Rasaiah  and  Stell.4 


5 


Using  quantum-mechanical  perturbation  theory  Bell  has  shown  that  the 
non-additive  three-body  potential  due  to  fluctuating  dipoles  and  multipoles 
may  be  expressed  as  a  function  of  the  intermol ecul ar  distances  and  the 
interior  angles  of  the  triangle  which  they  form 


< _ i _ . _ — - - — - i - - 

Rasaiah  and  Stell  recognized  that  the  corresponding  three-body  term  in  the 
perturbation  theory  of  polar  fluids  had  the  same  form 


U123  =  zil*2*3  W£]il2il3  (r12,r13,r23,cil’a2,:i3) 
where  Z£l£2£3  is  a  coefficient  determined  by  and 

Bell  derived  the  geometrical  factor  W  (r10,r, 

dipoles,  quadrupoles  and  all  triplet  combinations  of  them, 
his  work  to  include  octoooles  and  hexadecapol es ,  we  find 
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W144  =  105/16384  r^r ^3r 23  {105[9cos5a1+73cos3a1-12cosa1] 
+1080cos(a2-a3)[6-13cos2a-1-7cos4a1] 

+1188cos2(a2-a3)  [1 5cos3a-j  -COSa-j  ] 

+10296cos3(a2-a3)[l-3cos2ai ]+3861 0cosaicos4(a2-a3)} 

W 224  =  105/4096  ri2rl 3r23  1600+3267cos6a3-1728cos4a3+165cos2<*3 
+20cos(«1 -a2)[110cos»3-69cos3a3-297cos5a3] 

+70cos2  (a-|  -a2)  [33cos4a3+60cos2a3+35]l 

?44  334  344  444 

The  geometric  terms  W  ,  W  ,  W  ,  and  W  have  yet  to  be  obtained. 
Work  on  this  is  proceeding. 

£  £  £ 

To  determine  the  specific  coefficient  Z  1  2  3  in  the  three-body 
potential  we  made  use  of  the  spherical  harmonic  expansion  coefficients  and 
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the  method  outlined  by  Rasaiah  and  Stell  to  obtain 


113 

4  2 

=  pV/63, 

Z133  =  u2fi4/147. 

223  4  2 

=  ©V/175 

233 

4  2 

=  S7V/245, 

Z333  =  n6/343 

114 

4  2 

=  p  <t>  /81 , 

Z144  =  p2$4/243. 

7224  _  4  2  /9or 

Z  =  0  p  /225 

p  D  £ 

From  these  sequences  one  may  predict  the  values  of  L  1  23  even  when  the 
geometric  terms  beyond  W^  have  not  yet  been  obtained.  We  expect  that 
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To  complete  this  portion  of  the  work  we  needed  to  evaluate  the  remaining 
four  geometric  factors  to  obtain  the  effective  angularly-averaged  three-body 
potentials  which  contribute  to  the  term  of  0( a  )  for  all  binary  mixtures 
with  multipole  moments  through  the  hexadecapole  moment.  The  integrals  of 
the  form 

/g  i  £(  i^u-j  r]  2  ,UJ-j  dw^dr^  2 
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were  then  numerically  evaluated  using  a  method  previously  developed  by 
Rasaiah,  Larsen  and  S tell . ^  The  result  is  the  first  and  second  terms  for 
the  lambda  expansion  from  which  a  Pade  approximant  for  the  free  energy 
can  be  formed. 

In  order  to  extend  these  perturbation  methods  to  systems  with  non- 
spherical  cores  we  must  first  properly  characterize  the  reference  system. 

To  that  end,  we  proceeded  with  a  computer  simulation  of  linear  non- 
spherical  molecules.  The  objective  was  to  compute  an  ensemble  averaged 
angularly  dependent  distribution  function  in  terms  of  the  coefficients  in  a 
spherical  harmonic  expansion  of  the  distribution  function 
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From  these  coefficients  one  may  obtain,  by  the  orthogonality  of  the  spher¬ 
ical  harmonics  in  the  expansions  of  the  distribution  and  potential  func¬ 
tions,  a  very  simple  expression  for  the  term  of  0('O  for  a  linear  molecule 
with  a  quadrupol e-moment 
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where  is  the  diameter  of  one  of  the  atoms  in  the  diatomic  molecule  a  no  . 
is  the  reduced  density. 

Mor.te  Carlo  simulations  of  the  equilibrium  properties  of  a  nomonuclear 
diatomic  molecule  with  molecular  dimensions  characteristic  of  Cl-  were 
perfc  ~,iec  to  generate  an  accurate  set  of  coefficients  g  'r,  apoearing 

i  >.  '  J » i 
1  il 

if,  the  spnerical  harmonic  expansion  of  tne  radial  distr irtuicr.  'unction 
g'r,  Nearly  a  million  configurations  were  generated  in  calculating 

the  results  depicted  in  Fig.  1.  The  monvnts. 


wmcr  appear  in  the  term  of  0 ( • )  in  the  free  energy  decreased  stead’lv  by 

only  2. -a  over  the  last  BOOK  configurations.  These  calculations  are 

proper  :v  -ore  accurate  than  earlier  detensi nations^ of  g  (r),  even 

,  1  <  2r. 

thougr,  tne  programs  used  were  essential  1>  tne  same.  Tnis  is  primarily 
because  a  larger  number  of  configurations,  with  more  frecert  sam; linos  at 
let-  intervals,  were  generated.  !r  Tar’e  I,  the  newer  estimates  o'  tne 
ten,  of  2",  ■ )  are  compared  with  older  dare,  for  reduced  quadrupc i v  moments 
corresponding  the  critical  point  and  melting  point  of  Cl-  at  a  reduced 
density  of  0.5. 

Table  I 

Comparison  of  the  term  of  0(x)  obtained  by  numerical  integration 
or  (7)  using  the  data  for  g?-0{r),  g-?.(r),  and  g„„-(r)  from  the 
independent  studies  of  Street  and  T i f  des 1 ey  (S&T;  ahd  Eggebrecht 
and  Rasaiah  (E&R) 


Simpson 1 

$  rule 

Trapezoidal 

rule 

S&T 

E&R 

E&R 

S&T 

.812 

.283 

.328 

.373 

.414 

1  .4, 

.493 

.571 

.652 

.723 

1.975 

.686 

.798 

.907 

1.006 
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A  number  of  difficulties  became  apparent  from  the  discrepancies  in  the  table. 
The  sources  of  error  in  this  calculation  were  found  to  be:  1)  the  large 
fluctuations  in  the  moments  of  g^Cr);  2)  the  sensitivity  of  the  function 
to  the  method  of  integration.  The  greater  similarities  in  our  estimates  of 
the  integral  may  be  attributed  to  the  smaller  increments  used  in  our  version 
for  the  radial  distance. 

II.  Adsorption  of  Dipoles  and  Ions  at  a  Wall  in  the  Presence  of  an  Electric 
Fi  eld 

While  our  understanding  of  the  behavior  of  liquids  has  improved 
dramatically  over  the  past  ten  years,  the  physics  and  chemistry  of  surface 
phenomena  is  still  at  a  primitive  stage.  One  of  the  outstanding  problems 
in  this  area  is  the  elucidation  of  the  density  profiles  of  ions  and  dipoles 
at  a  charged  surface;  a  detailed  knowledge  of  this  is  fundamental  to  our 
understanding  of  phenomena  ranging  from  the  behavior  of  electrodes  in 
batteries  to  selective  adsorption  and  transport  across  membranes.  It  is 
convenient  to  study  ions  and  dipoles  at  a  surface  separately  before  trying 
to  gain  insight  into  the  behavior  of  a  mixture  of  these  two.  We  have, 
therefore,  embarked  upon  a  program  in  which  the  following  problems  will  be 
studies  in  sequence. 

1)  Dipoles  at  a  wall  in  the  presence  of  an  electric  field. 

2)  Ions  at  a  charged  interface. 

3)  Dipoles  and  charges  at  a  wall  in  the  presence  of  an  electric  field. 

Investigation  of  the  first  problem  has  progressed  very  rapidly  in 

collaboration  with  Dr.  Dennis  Isbister  and  John  Eggebrecht.  Professor  Noel 
Thompson,  a  visitor  on  sabbatical  leave  from  Australia,  Dr.  Dennis  Isbister, 
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and  I  have  made  progress  on  the  second  problem,  while  work  on  mixtures  of 
ions  and  dipoles  at  a  wall  has  only  just  begun.  In  what  follows,  we  provide 
a  brief  description  of  our  work  on  dipoles  at  a  wall. 

The  effect  of  a  wall  is  obtained  by  allowing  a  mixed  system  of  at  least 
two  species  1  and  2  to  change  in  such  a  way  that  the  density  of  one  species 
(2  say)  goes  to  zero  accompanied  by  an  increased  in  diameter  R^.  It  is  the 
growth  of  this  single  particle  (2)  that  emerges  finally  as  a  wall,  and  if 
the  wal 1 -particl e  had  a  dipole-moment,  the  zero  density  -  infinite  radius 
limit  could  be  made  to  produce  an  electric  field  as  well  emerging  at  any 
chosen  angle  from  the  wall.  The  density  profile  of  each  fluid  species 
near  the  wall  is  given  by 

c,(r)  =  lim  lim  p °[h21 (p ,^2)+l ] 

where  c^°is  the  bulk  density  of  species  1,  and  h^fr)  is  total  wall-particle 
correlation  function  which  has  the  usual  invariant  expansion. 

Our  initial  investigation  of  dipoles  at  a  wall  used  the  linearized 
hypernetted  chain  (LHNC)  closure  for  the  wal 1-particle  and  particle-particle 
interactions . 

It  is  known  that  this  approximation  leads  to  fairly  accurately  bulk 
properties  for  dipoles.  We  found  that  the  contact  values  for  the  wall- 
particle  correlation  function  are  much  higher  than  those  obtained  in  the 
mean  spherical  (MS)  approximation/ 

The  linearized  hypernetted  chain  closure  is  obtained  by  taking  the 
first  term  in  the  expansion  of  the  angular  part  of  lnth^+l)  which  appears 
in  the  hypernetted  chain  (HNC)  approximation  for  the  direct  correlation 
function  .  If  the  second  term  is  retained  as  well,  a  more  accurate 


12 


description  of  the  system  results.  This  approximation,  known  as  the 

8  9 

quadratic  hypernetted  chain,  QHNC,  ’  reveals  a  new  phenomenon  that  is 
absent  in  the  MS  and  LHNC  approximations.  This  is  the  change  in  the 
density  of  the  fluid  in  an  open  system  when  the  electric  field  is  turned 
on,  and  is  called  el ectrostriction .  The  leading  term  in  the  relative 
change  in  density  Kh  at  an  infinite  distance  from  the  wall  was  found,  in 
the  QHNC  approximation,  to  be  given  by  a  term  of  0(E  ) 


k  (2)  _  Ap_ 
h  PjU 


u-i) 


2  E 


2 


Q 


where  E  is  the  electric  field,  Q  the  inverse  compressibility,  c  the  di- 

4-  2  o 

electric  constant,  and  y  =  y-  m^  6,  with  p  =  1/kT.  Approximations 
beyond  the  QHNC  theory  systematically  generate  term  of  higher  order,  so 
that  this  approximation  must  contain  the  complete  electrostriction  term 
of  0(E  )  in  the  HNC  approximation.  However,  a  discrepancy  in  was 
found  on  comparing  our  result  with  the  thermodynami cs  of  electrostriction 
in  an  open  system.^  Since  the  HNC  approximation  and  all  of  the  theories 
derived  from  it  ignore  bridge  diagrams,  the  discrepancy  in  the  term  of 
0(E  )  between  the  QHNC  result  and  thermodynamics  must  come  from  neglect  of 
these  diagrams.  The  lowest  order  bridge  diagram  was  then  calculated  and 
found  to  contribute  significantly  to  electrostriction.  Adding  this  to  the 
QHNC  result  of  0(E  ),  we  obtained  (resolving  the  above  discrepancy) 


K  b(c-1)2  5(t-l)26,yl  E2 

h  ^24-rtp^uy  1  28np  o  jq 

which  shows  that  the  bridge  diagrams  cannot  be  ignored  in  calculations  of 
the  local  density  at  an  infinite  distance  away  from  the  wall.  It  appears 
likely  that  they  affect  the  density  profiles  also  when  the  dipoles  are 
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closer  to  the  wall.  This  is  an  important  result,  because  all  of  the 
theoretical  studies  of  the  ordering  of  dipoles  and  ions  at  a  wall  that 
have  been  carried  out  so  far  ignore  these  bridge  diagrams. 

Our  study  of  the  LHNC  approximation  for  the  wall -particl e  distri¬ 
bution  function  leads  to  the  constitutive  relation 

?(“>)  =  (e-1)  ?(®)/4u 

between  the  polarization  denisty  ?(<*=)  and  the  electric  field  f  with  the 
same  dielectric  constant  c  as  that  obtained  in  relating  e  to  the  two- 
particle  correlation  function  of  a  dipolar  fluid  in  the  absence  of  an 
electric  field.  The  QHNC  approximation,  however,  yields  in  addition 
non-linear  terms  in  the  electric  field  for  the  polarization  density. 

We  plan  to  extend  this  work  to  mixtures  of  dipoles  in  order  to  study  the 
selective  adsorption  of  dipoles  near  an  electrified  wall. 
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Fig- 1  The  coefficients  of  the  spherical  harmonic  expansion  of  the  angle  dependant 
distribution  function  after  two  million  configurations. 
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A  theory  of  elect  restriction  which  follow  v  from  studies  of  dipolar  ordering  at  an  electrified  wall  is  discussed  in  the  qua¬ 
dratic  hypernetted  chain  approximation,  bridge  diagrams  for  the  wall  particle  correlation  functions  contribute  significant¬ 
ly  to  electrostriction  even  to  lowest  order  in  the  electric  field.  The  form  ol  the  constitutive  relation  between  the  polariza¬ 
tion  density  and  the  field  in  strong  fields  is  discussed. 


I .  Introduction 

The  ordering  of  dipoles  at  a  wall  from  which  an 
electric  field  emerges  has  been  studied  in  the  mean 
spherical  (MS)  (1  ]  and  linearized  hypernetted  chain 
(LHNC)  1 2 )  approximations.  Wc  icport  here  that  elec¬ 
trostriction  appears  as  an  added  featuic  when  the  the¬ 
ory  is  carried  beyond  the  LHNC  approximation  for 
the  wall  particle  correlation  functions,  leading  to  a 
molecular  theory  of  this  phenomenon  in  w  hich  graphi¬ 
cal  analysis  and  integral  equation  approximations  that 
are  ubiquitous  in  the  theory  ol  fluids  may  be  exploit 
ed.  Of  the  hierarchy  of  approximations  generated  by 
the  hypernetted  chain  equation  (UNO.  we  find  that 
the  quadratic  hypernetted  chain  (01  INC)  approxima¬ 
tion  is  the  first  to  predict  electrostriction.  The  leading 
term  in  the  relative  change  in  density  an  infinite  dis¬ 
tance  away  from  a  flat  wall  is  of  0(/  *  I.  where  E  is  the 
magnitude  of  the  local  field.  Approximations  beyond 
the  QHNC  theory  systematically  generate  terms  of 
higher  order,  so  that  the  QHNC  approximation  also 
contains  within  it  the  complete  electrostriction  term 
of  0( E*)  in  the  HNC  approximation. 


All  of  the  theories  of  electrostriction  which  have 
i hoi r  genesis  in  the  HNi  .  ,-proximation  ignore  the 
bridge  diagrams,  but  w  e  find,  by  comparing  the  QHNC 
theory  with  the  thermodynamics  of  electrostriction 
|.M.  that  these  diagrams  must  contribute  significantly 
even  to  lowest  order  in  /  This  leads  us  to  believe  that 
they  play  an  equally  important  role  in  determining  the 
local  density  when  the  dipoles  are  closer  to  the  wall. 

We  lepoit  hete  the  baicst  outline  of  our  molecular 
1  booty  lot  an  open  system  |4]  in  which  contributions 
Horn  the  QIINC  approximation  to  0(f.  )  and  from 
the  bridge  diagrams  to  the  same  order  in  the  electric 
held  and  to  lowest  older  in  the  density  and  the  dipole 
moment  aie  evaluated  analytically.  Our  theory  can  al¬ 
so  be  extended  to  higher  order  in  /:  but  this  will  not 
be  pursued  hete.  The  thermodynamic  discussion  of 
Kirkwood  and  Oppenheim  |3 )  however  is  limited  to 
deriving  the  elect  tost  net  ion  term  of  0(/f*),  since  it  is 
based  on  a  thermodynamic  perturbation  theory  of 
first  order  in  /  2 .  that  uses  as  reference-state  input  the 
linear  constitutive  relation  between  the  polarization 
density  Pi'*')  and  the  local  electric  field  E  for  a  field- 
independent  e. 


Volume  79.  number  2 


CHFMK'AI  PHYSICS  LETTERS 


15  April  1981 


Pi00)  =  U  1 )  F(°°)/4n.  (  1  I 

Turthei  details  of  our  work  and  extensions  to  it 
will  be  presented  elsewhere  |4|  but  we  remark  tluu 
oui  studs  yields  the  constitutive  relation  with  the 
same  value  of  e  as  that  obtained  from  from  relating  < 
lo  the  two-particle  orientational  correlation  function 
in  the  absence  of  a  field,  and  also  provides  a  basis  toi 
determining  non-linear  terms  in  the  electric  field  loi 
the  polarization  density. 

2.  Electrostriction  in  an  open  system 

We  employ  the  technique  discussed  in  refs.  1 1 .2 1  to 
generate  an  electric  field  by  taking  a  mixture  of  dipo¬ 
lar  hard  spheres  to  the  zero-density  limit  of  species  2 
with  an  attendant  increase  in  the  radius  A’i  to  infinity . 
under  the  constraint  that  its  dipole  moment  divid¬ 
ed  by  the  cube  of  the  radius  of  the  excluded  volume 
A’,,  =A-,  +A,  is  a  constant  A'p.  The  electric  field  /  s 
emerging  from  the  wall  is  related  to  by  |1 .2) 

Ft  =  /■(>(.*  cos-Oj  +  1 1*  2  <*t.  1 2  l 

wheiet't  is  a  unit  vector  dependent  upon  0-,  |l.2|. 
the  angle  which  the  wall  dipole  makes  with  the  wall 
normal  n  Note  that  Ft  is  constant  when  the  wall  di 
pole  orientation  6 1  is  fixed.  The  relationship  between 
Ft  and  the  Maxwell  field  F  in  the  fluid  is  of  the  toim 

F  =  |3  (2e  +  I  )J  Ft  +  terms  non-linear  in  Ft.  (7i 

ileie  the  non-linear  terms  aie  of  magnitude  const  /  2 
+  .  when  Oj  =0,and  c  is  the  dielectric  constant  ol 

the  Hind  at  zero  field.  In  the  notation  ot  ref.  |2  | .  the 
wall  particle  correlation  function  has  the  invan.im 
expansion 

h  t|  1  - .  A  y . H  \ )  =  h 2j  (r )  +  /r  t]  lei  /2( 2 . 1 ) 

+  ^(2.  1)  ♦ _  (4) 

where  A (2, 1  )  =Jt  •  jj  and  77(2, 1 1  =  St-(3n«  U  t  ■  s, 
in  which  .<->  and 5 1  arc  unit  vectors  in  the  directions  ot 
the  wall  dipole  and  fluid  dipole  respectively  ,r  is  the 
distance  of  a  fluid  dipole  (of  orientation  ft ( 1  from  the 
wall,  and  U  is  the  unit  tensor.  In  this  wall  limit .  A-. 

•►oo_ 


where  /t^lc)  is  short -tanged  and  A  is  a  constant  re¬ 
lated  to  the  elect t ic  field  F>  (see  section  7).  Hence 

lim  j  (cl  =  .'A  (M 

The  oiigin  ot  eleetiostnctioii  in  our  theo¬ 
ry  is  the  coupling  that  can  exist  between  /i(J  (r)  and 
hi]  (c):  this  feedback  begins  w it  1  >  the  01  INC  appiox- 
utiation 

We  define  the  electros!  net  ion  ettect  A;,  as  the  tela- 
tive  change  in  density  ot  the  bulk  thud  when  the  dec- 
trie  field  is  sw  itched  on.  so  that 

K/,  =■  Ap/pf|  =  li\ |t«’.  Fs ),  ( 2) 

where 

/i$1u.f2)  =  s2  ‘J/i^t.-.Fs.n,  win,  <m 

il  ~  4rr  for  dipolat  hard  sphetes  and  p"  is  the  density 
of  the  bulk  fluid  when  /  -  -  0.  Il  A  is  the  correspond¬ 
ing  asymptotic  limit  of  the  angulatly  avetaged  diteci 
cort  elation  I  unction  <s|  l :.  /  .  .12  (  I.  we  find  |4|  that . 
tn  the  absence  of  molcenlai  poiat  tz.i1' tit t x 

F/,  =  A,  Q. 

where  (7  is  the  invetse  comptes'ibibty  "•  the  tluid. 

(7=1  ( Pj‘  -S2“  >  /  <  il.)  dfl|  dft,  dz  (101 

and  t'ult'.  12,  .  ftr  I  is  the  dnee!  coiie'atton  function 
of  the  bulk  fluid.  The  thermndy  natm.  theory  o!  elec¬ 
tros!  fiction  |7|  tot  an  open  systetti  gives  A/(  toOiA ’  I  as 

A/,2'  =  (  (3'Htt  )  ( <‘*e  V,'i/  2  ('.  (Ill 

where  (3  =  1  Ik T.  k  is  Holt/mannN  constant .  7  is  the 
absolute  temperature  and  the  stipct script  (2)  means 
the  let m  of  0 (/  2  ). 

3.  The  OHNC  approximation  for  electrostriction 

By  considering  the  asymptotic  limit  of  (  ^(c)  in  the 
invariant  expansion  tot  the  wall  panicle  direct  corre¬ 
lation  I  unction  (i|  (c.  I .  s  12  j  two  find  |4|  that  A->j  is 
related  to  /  s  in  the  <>ll\(  apptoximat toil  by 


O  t  „i  -  t.P  t  - 
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,/:/,( 3  cos 202  +  1)1'2  =  A'21  l2C>+(2A:llpV/ifii  > 

+  t>  (  A,^?,)  3A,,/(1+A,,)1.  (12) 

where  the  (3,(  )  functions  are  defined  itr  ref.  |2|. 

and  the  constant  A  is  related  quite  generally  (i.e.  not 
just  m  the  Q1INC  theory)  to  the  fluid  dipole  moment 
in,  by  |4  6 1 

3t  =C>+(2A'up,i/?f1)  V  A- ■A.'uP1, ’«•},).  (131 

whole  r  =  47r«t?p?0/9 .  The  dielectric  constant  e  of 
the  fluid  is  also  quite  generally  given  by  |6) 

e  =  Aup?A-},).  (141 

We  also  find  from  the  z  -*  00  limit  of  the  QHN(  ap¬ 
proximation  for  oji  £/  .Oj )  that 

A(l2’  =4  A 5i  (3  cost 0 2  +  1)  (15) 

and  on  using  (2),  (3),  (9)  and  (12)  ( 14)  we  have 

A,1,2’  =  ((3'24rrpO».)(e  -  l)2  E2/Q.  (Us) 

where  the  superscript  (2)  means  again  the  term  of 
O (/  -).  Note  that  in  (16)  A.'/2’  is  independent  of  the 
inclination  0-,  of  the  wall  dipole.  Consistency  with  the 
thermodynamic  theory  (11)  demands  that 

pV  Oe  Vi  =(f  D2/3y.  (IT) 

This  relation  is  indeed  satisfied  by  the  simple  Debye 
equation  (e  -  1  )He  +  2)  =  r.  but  the  more  exact  ex¬ 
pression  [7.8]  (for  arbitrary  y  and  p* 1  -*0) 

tc  1 )  ( e  +  2 )  =  v  g  .>,3+-  <<*> 

introduces  a  discrepancy  of  0(>,2)in(17).  Since  lnglier 
order  approximations  derived  from  the  HNC  equation 
do  not  produce  additional  contributions  to  A'/,  of 
(Ml  -  ),  the  apparent  inconsistency  lies  in  ignoring  the 
bridge  diagrams. 


4.  Bridge  diagrams  of  0(£2) 


The  bridge  diagram  of  OlA2)  and  to  lowest  order 
in  the  fluid  density  can  be  represented  graphically , 
when  r  —  by 


B*(2. 1 )  = 


(l‘>) 


where  c - -  •  =  3A'„  D(Z.  i).  and 

■>  “ 


is  the  three-particle  direct  correlation  function  r3(A) , 
,V3 ,  A’4 )  for  the  fluid  with  ,V,  2  (r,,  Sll  1.  li‘ (2. 1 )  con¬ 
tributes  to  the  wall  -  particle  direct  correlation  func¬ 
tion.  In  (19)  the  open  circle  2  is  the  wall  dipole  root 
point,  the  half  black  circle  »  signifies  angular  integra¬ 
tion  over  the  orientations  of  the  bulk  fluid  dipole  1 
and  the  black  circles  (field  points)  lepresent  spatial 
and  angular  integrations  of  fluid  dipoles.  To  evaluate 
(19)  to  lowest  order  in  P]'  and  nq  we  use  the  low-den¬ 
sity  limit  in  terms  of  Maver/ functions 

c3a, .  *3.  -V4  )  ^  )\X  j .  -V3  )  /(*3,.V4  )  f(X4.Xx ) 

(20) 

and  retain  only  the  n  =  1  term  in  the  perturbation  ex¬ 
pansion  [8.9) 

J\Xj.  Xj)  =/f|(r,;  ‘ 

+  [1  +/n(r„)l  2  [0m]DUJ)irl\H /t\\  (21) 

n  - 1 

for  the  Mayer  /  function,  where  fn  (rJ(  I  is  the  corre¬ 
sponding  Mayer  function  for  the  reference  hard-sphere 
system.  The  bridge  diagram  is  now  evaluated  analyti¬ 
cally  using  ilankel  transforms  and  we  find  eventually 
that  when  c  -*  00 . 

Ha),(oo)-y  rfs(e  1  )2  fifA 2  np(,’ .  (22) 

in  which  the  contribution  from  the  three -particle  cor¬ 
relation  function  appears  as  the  third  virial  coefficient 
for  hard  spheres!  Adding  this  to  ( 16).  alter  using  (9), 
we  have 

*(2)J-J3(e  1>2  5(e  1  )20v 

h  L  24irp(>  v  1 28trp9 

where  —  emphasizes  that  our  calculation  of  the  bridge 
diagram  of  0(/'2)  is  correct  only  to  lowest  order  in 
p^  and  m  j .  This  is  sufficient  however  to  test  the  mo¬ 
lecular  theory  against  thermodynamics  using  (18)  for 
the  dielectric  constant.  Instead  of  ( 1 7)  we  now  need 

p?3e/dp?~(e  l)2 '3v-  £(6  1  l2r.  (24) 

which  is  consistent  with  1 18)  to  0(  r2).  We  have  thus 
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shown  that  the  bridge  diagrams  make  significant  con¬ 
tributions  to  the  density  when  and  it  appears 

likely  that  they  affect  ilie  density  profiles  also  when 
the  dipoles  aie  closer  to  the  wall.  Theoretical  studies 
of  the  ordering  of  dipoles  and  ions  at  a  wall  that  have 
appeared  so  lar  ignore  these  bridge  diagrams. 


5.  The  polarization  density  P{°°,  A'i )  in  an  open 
system 


From  the  definition  ol  Pyz.  F\ )  in  the  grand  ensem¬ 
ble  |t>)  we  have 

P(z.E->  1  =  («tjPj>  12 1  (25) 

In  our  discussion  Fi  is  independent  of  z  [see  (2)| . 
Using  the  invariant  expansion  for  /tjjltt.Ifj.  fi, ).  we 
find,  after  doing  the  necessary  angular  integrations, 
that 

P(°°.  A\ )  =  Wjp'i'A  >|t.'  cos-0,  +  H1  2  e,.  (2o) 

which  is  independent  of  any  approximation  for  the 
wall  particle  correlation  function!  Employing  the 
QHNC  approximation  lor  K^\  given  in  (12), 

P^.t.  :  i  =  (3p(,'miF;|2c>+(2>v  „P?A'i,) 

+  (>  (  A  , ,  pVA’ii  >  -’A/,/  (I  +  K/,  )|  1  (271 


=  |(e 


X 


1  >  4tt)  l-'/r  ,  (2f  +  1)1 

Kh  €  l'|  1 

.1(1  +  Kh\  2e  +  1J  ' 


(28) 


When  K)t  =  0  (i.e.  m  the  absence  of  elcctrostriction  as 
in  the  MS  and  L.HNC  approximation)  and  when  all  the 
higher  coefficients  of  /:?  (rr  >  3)  in  (3)  arc  assumed 
zero,  we  recover  the  constitutive  relation  ( 1 ).  Con¬ 
versely  when  K tl  and  the  higher  coefficients  in  (3)  are 
not  zero,  non-linear  terms  in  the  polarization  density 
will  appear.  Substitution  of  ( 1 5)  and  (9)  in  (28)  follow¬ 
ed  by  expansion  of  the  denominator  draws  out  the 
term  of  0(/  s  )  in  the  QHNC  approximation  for  the  po¬ 
larization  density  which  is 

PV\\\(  {oc ■  F  ]  =  (<c  '  1  )/4ir|(3f2/(2e  +  l)| 


X 


e  1 
2e  +  I 


(29) 


192 


If  the  non-linear  terms  in  (3)  are  neglected 
/V>IINC<~  /  l(f  l)  4rr|/ 


1  + 


24~p')li' 


(e  1)? 

2f  +  1 


/- 

L> 


+  0(/  5i. 


(30 1 


A  more  complete  calculation  of  this  term  in  the 
QHNC  approximation  would  teipnie  the  deteimina- 
Hoi.  ol  i he  coefficient  of  the  /  term  in  ( 5 1  when  6  . 

“  0.  An  exact  calculation  of  the  term  ol  OU  '  1  w  ould 
tevjune  m  addition  that  the  contributions  ot  the  bridge 
diagiams  to  elcctiostriction  Jsee  i22l|  and  to  the  rela¬ 
tionship  between  K  s,  and  /  ■,  he  included. 
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The  adsorption  of  dipoles  at  a  wall  in  the  presence  of  an 
electric  field:  The  RLHNC  approximation50 
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The  adsorption  of  dipolar  hard  spheres  in  the  presence  of  an  external  electric  field  has  previously  been 
studied  within  the  context  of  the  mean  spherical  approximation.  In  order  to  quantify  the  significance  of  the 
physical  trends  found  above,  the  problem  is  solved  within  the  higher  order  closure  rules  afforded  by  the 
linearized  hypernetted  chain  approximation.  I  xpressions  for  the  reduced  dipole  moment  and  the 
electric  field  strength  are  derived  using  only  the  asymptotic  forms  of  the  direct  correlation  functions  It  is 
found  that  the  favorable  orientational  correlations  between  the  dipolar  hard  spheres  and  the  wall  are 
underestimated  by  the  mean  spherical  approximation.  This  is  emphasized  in  the  enhanced  adsorption  of  the 
dipolar  species  (at  the  wall  itself)  for  dipoles  oriented  close  to  the  direction  of  the  field.  However,  the 
nonphysical  features  of  the  mean  spherical  approximation  (manifested  in  the  negativity  of  the  density  profile) 
are  not  fully  rectified  by  the  use  of  the  linearized  hypernetted  chain  approximation. 


I.  INTRODUCTION 

The  adsorption  of  dipolar  molecules  at  a  wall,  in  the 
presence  of  an  electric  field,  is  of  interest  in  the  study 
of  electrode  and  membrane  phenomena.  Here  the  ad¬ 
sorption  phenomenon  is  delineated  by  the  distribution  of 
molecules  at  a  particular  orientation  12,  and  distance  v 
from  a  hard  planar  wall  p,(z,  E2,  12,).  The  electric  field 
E2  emanates  from  this  wall,  the  declination  of  the  field 
with  respect  to  the  wall  being  allowed  by  its  nonconduc- 
tive  properties.  Isbister  and  Freasier1  have  investi¬ 
gated  this  problem  for  hard  dipolar  spheres  against  a 
hard  wall  using  the  mean  spherical  approximation 
(MSA).  Their  results  for  the  density  profile  p,U,  Eg,  12,) 
of  dipoles  are  of  great  interest  even  though  they  suffer 
from  the  defect  that  the  wall  particle  density  profile 
p,(z.  E2,  12,)  assumes  negative  values  at  certain  relative 
orientations  of  electric  field  (E2)  and  dipole  moment  m, 
of  the  particles  in  the  fluid  (the  dipole  orientation  is  de¬ 
noted  here  by  12,).  However,  the  argument  leading  to 
the  electric  field  at  the  wall  is  not  swayed  by  the  approx¬ 
imation  used,  and  may  be  employed  with  more  accurate 
theories  such  as,  for  example,  the  linearized  hypernetted 
chain  (LHNC)  approximation.  While  these  theories  could 
be  expected  to  produce  better  results,  they  do  suffer  from 
the  necessity  of  employing  numerical  methods  to  a 
greater  extent  than  is  needed  to  determine  the  density 
profiles  in  the  mean  spherical  approximation. 

This  paper  is  devoted  to  a  study  of  the  wall-particle 
density  profile  using  the  linearized  hypernetted  chain 
approximation2  for  the  wall  particle  and  particle-parti¬ 
cle  interactions,  except  that  the  effects  that  are  inde¬ 
pendent  of  the  orientations  of  the  electric  field  and  the 
fluid  dipole  are  treated  exactly.  In  our  study,  these 
are  the  interactions  between  the  hard  cores  in  the 
fluid,  and  also  the  interactions  between  these  cores  and 
the  hard  wall.  We  call  this  the  renormalized  linearized 


*’ Extracted  in  part  from  the  M.S.  (Chemistry)  thesis  of  .1 
Eggehrecht,  t  niversitv  of  Maine  (1980). 

"Present  address:  Department  of  Chemistry,  Faculty  of 
Military  Studies,  I'niversity  of  New  South  Wales,  RMC, 
Hunt roon,  ACT  2600,  Australia. 


hypernetted  chain  approximation  (RLHNC)  after  the 
nomenclature  introduced  by  Stoll  .md  Weis.3  Our  re¬ 
sults  for  this  theory  are  an  improvement  over  the  cor¬ 
responding  MS  approximation  when  the  system  is  char¬ 
acterized  by  (a)  a  reduced  fluiddensity  p*  -p,/? ,,of  0.  5711. 
(b)  a  reduced  dipole  moment  m*  m,  v  fcTRJ,  of  v  0.5 
(or  equivalently  a  reduced  temperature  T*  ^  1  in'2  of  2), 
and  (c)  a  reduced  external  electric  field  E *  =  EZR\ i  m, 
of  8/3.  In  contrast  to  the  MS  approximation,  the 
RLHNC  wall-particle  density  tunctions  p,(z,  E2,  f2,)  are 
only  marginally  negative  near  the  wall  for  a  dipole  ori¬ 
entation  in  direct  opposition  to  the  field  (see  Fig.  2). 
These  functions,  however,  can  become  negative  over  a 
larger  distance  z  from  the  wall  when  the  reduced  dipole 
moment  is  increased  to  1.0  without  altering  the  reduced 
electric  field  or  reduced  density  (f  Jg.  5).  This  sug¬ 
gests  even  higher  order  terms,  beyond  the  RLHNC  ap¬ 
proximation,  must  be  included  in  (he  theory  when  the 
dipole  moment  ///,  and  the  external  electric  field  Ez  are 
both  large. 

II.  GENERAL  THEORY 

The  technique  of  producing  a  wall  next  to  a  fluid  bj 
taking  the  limiting  behavior  of  a  binary  mixture  (with 
densities  p,,  p2  and  radii  R,,  Rz)  detailed  by 

lim  lim  (2.1) 

^2*  •  02~  0 

is  well  known.4  Isbister  and  Freasier1  have  extended 
this  to  introduce  an  electric  field,  as  well,  by  consid¬ 
ering  the  corresponding  limit  for  a  dipolar  mixture  un¬ 
der  the  restriction  that  the  dipole  moment  w2  of  par¬ 
ticle  2,  which  eventually  becomes  the  wall,  divided  by 
the  cube  of  the  radius  of  the  excluded  volume  R2,  =R2 - R, 
is  a  constant  (hereafter  called  E0): 

lim  i»i/R21  =  £0  .  (2.2) 

8  2"  “ 

In  taking  these  limits,  in  the  specific  order  pt-  0,  R2 
— 00 ,  the  volume  of  the  system  is  allowed  to  grow  faster 
than  R2,  keeping  p,  constant  through  the  constraint  that 
PzR\~  0. 

The  magnitude  and  direction  of  the  electric  field  E2 
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llu-  ‘‘wall  dipole.  -  which  lias  receded  to  a  distance  -  ' 
from  the  wall.  At  that  distance,  the  vector  r£,  become? 
perpendicular  to  the  wall,  and  the  electric  held  which 
makes  an  angle  o,  with  respect  to  this  normal  is  riven 
by  the  appropriate  solution  to 


coso,  r,  •  »■,,  - 


2  cose 


(3  cos  t,  .  1 ) 


• 


(2.  101 


l  lu.  ’.  1  he  eiev'.n  held  ai.glv  ai  Ira  wall  nl  itied  against 

*.1k  uwlinaiiun  t>  >1  ‘Vail  ii;  ale”  c-  '  resides  at  minus 
inlim:  and  :'r,xj..ccx  ‘he  electric  :,eld  £  .  IVa  angle  S,  is  the 
inclination  at  a  m;.ole  in  tile  i .  mud. 


follow,  when  this  limit  is  applied  to  the  dipole-dipole 
Interaction  energv 


Tlie  solution  to  this  equation  pices  the  direction  ol  the 
electric  field  uniquely,  in  terms  of  the  orientation  ol 
the  wall  dipole  at  minus  infinity  (see  Tip.  1). 

The  density  profile  E£.  11,1  of  dipolar  molecule? 

at  a  distance  .  trom  tin  wall  depends  on  the  dipole  it", - 
enlalion  C!lt  ■.«.  the  direction  ol  the  electric  tieln  E:  and 
on  other  parameters  of  the  system.  It  is  obtained  :  i‘<  ■  s : . 
the  relation1 

e2,  inn  *  am  ,t,  i;-.  <-•,)  •  1  f  .  12.  if 

W«  -  *•  |  i>2-  “  \ 

where  , > ,  is  the  bulk  density  ol  species  1  and 
h,l  is  the  tola!  correlation  (unction  ol  sprites  1  .,:w  2 
a  binary  mixture.  The  latter  is  the  solution  to  me  Urn- 
stein-  Zermke  relation 


,  1  V  , 

h2i  -  c2,  -  ^  4  (V,  . 

q  >=i 


L>(2.  1)  . 

(2.3) 

where  "  j  f/r3</Sl3  is  . 
angular  integrations. 

^7 21  < 7*21 »  ^2*  fl,)  -  (<2i(>'2j 

**  2  » 

(2.4) 

c21(r2,.  fl2.  11,)  c2,()  2, 

Here 


s,  and  s2  are  unit  rectors  in  the  directions  oi  the  dipoie 
moment  vectors  m,  and  m2.  respectively.  (3r2,  r2,  -  U) 
is  the  dipole-dipole  interaction  tensor,  in  which  ra  is 
a  unit  vector  in  the  direction  along  the  line  joining  par¬ 
ticles  1  and  2.  and  U  is  the  unit  tensor.  The  potential 
may  be  written  in  terms  ot  the  electric  field  Es,  pro¬ 
duced  at  the  location  of  particle  1  by  tin-  second  parti¬ 
cle  2. 


(2.  12! 


1.2.  13) 

),(r21)  M2.  1)  -  r£,<r21)/M2.  D  . 

(2. 14) 

and  using  Werthcim’s  multiplication  table"  in  Fourier 
space  lor  the  angular  integrations,  tlie  Ornstein-.  crniki 
relation  reduces  to  three  equations 


'‘’2i  -  r2i  '  Ix'V 4  * 


[2.  15' 


1 2i(r2i ,  1; ,,  1;  | )  -  —  )/,  i  s  i  •  e2  . 

w  he  re 

E2  -  -t*  e.  ■ 

>■  f,  - 

In  Eq.  (2.6' 

e2  ■  f 3 » ■  2 1  >"2i  -  U>  ■  x2  -  (3r,,  cose.  -  ? ,) 


(2.5) 


(2.6) 


(2.7) 


i<‘l i  -  4  • L  •  H-.  -  4  •  4  ■  4  •  .  (2.  io» 

°  >4 


*  4  -  4 . 4,t 


(2.  17) 


4i  -  cZ,  •  „  .  p,(2/i 

J  >»i 

when  *  !  i/r,  is  a  convolution  involving  only  spatial  in 

tegrations.  In  Eqs.  (2.  16)  and  (2.  171 


and  e2  is  the  angle  which  the  diiiole  embedded  in  particle 
2  makes  with  r21  (see  Fig.  1).  Since  the  magnitude  of 
e2  is  (3  cos2b2  -  1)'  2. 


iC(r)  =- 4;,(  1 )  -3  f  lb  X"1  4f( s)  12.  1H) 

and  its  inverse  is 

3  C 


—■ y2  (3  cos2b,  .  1 ) 1 '  :r2  . 


4„(r)  -4e(r)  -  ^  f  4„(s)  s*,is  . 


12.  19) 


(2.8) 


where  c2  is  a  unit  vector  in  the  direction  of  E2.  Changing 
variables  to  r21  -  Rz  -  z,  and  taking  the  wall  limit,  the 
electric  field 


with  similar  expressions  for  r£e(r)  and  cjg(>")  in  which 
a,  (j  =  1  or  2.  Inside  the  hard  cores 


4„(r)-5  • 


(2.  20) 


E2  : E0(3 cos2b2 -  l)WJr2  . 


(2.9) 


which  shows  that  E2  is  independent  of  the  distance  z  from 
the  wall,  but  that  its  magnitude  is  determined  by  the 
strength  (through  E0 )  and  the  direction  (through  02)  of 


*£«(»■)  =,'Sp(»>  -°  .  =  /(«.  1  • 

so  that  from  Eq.  (2.  18) 

*2.(r)--  3At"aS.  >■«„„,  (2.21) 

where 
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A„8  I  Ifals)  .  (2.22) 

•SjD 

Following  We rtheim. 3  the  relations  for  Jt£,(r)  and 
/ift(r)  may  be  uncoupled  by  taking  linear  combinations 

'<;„(>-)  "  ['<«>•>  *  “•  >  I  3Aap  .  (2.  23) 

h'atir)  -  fe(r>  -  *£„(>•>  I  ■  <2.  24) 

with  analogous  expressions  for  c*s( »•)  and  e"e(r).  The 
closure  conditions  (2.20)  when  applied  to  Eqs.  (2.23' 
and  (2.24)  are  equivalent  to 

/i*ai(r)--l  .  r  /ia,  .  (2.25) 

On  using  the  linear  combinations  h'atl(r)  in  Eqs.  (2.  16) 
and  (2.  17),  one  finds 

*;,(<•>  -  r*,,(r)  -  £  A'»  *5  '4  *  '■>>  •  <2-  2f>) 

r=l 

where  p  *  -  2p,,  p"  — Equations  (2.  15)  and  (2.26'  re¬ 
semble  the  Grnstein-Zermke  relations  for  a  binary 
fluid,  in  which  the  molecular  interactions  are  spherical¬ 
ly  symmetrical. 

In  the  limit  A2-  -i ,  Eq.  (2.  19)  together  with  Eqs. 
(2.23)  and  (2.24)  yields 

/i?,U)  =  [2/i;,(£)-//ilU)*3|A'!I  .  (2.27) 

As  p2- 0,  Eqs.  (2.  15)  and  (2.26)  reduce  to 

//!,(>•)=  <1i(r)  *  P,  *  <•;,  (2.28) 


(<2i(r)  -  +  b’tipt  ('21  *  (-ii  • 

which  can  be  written  in  bipolar  coordinates  as 


mined  by  the  interactions  between  the  particles  in  the 
bulk  fluid.  An  analogous  equation  can  be  written  for  the 
angularly  independent  part  of  the  wall-particle  total  cor¬ 
relation  function 

h'aU)  -  <‘21(c)  •  2*p,|z  (B'W  -B‘U)|-  [D‘(*)-U'U)\ 

■  Pdv>i!1(  v)lR5K)-Bs(U-  v|)||  ,  (2.35) 

where  /!'(.  '  and  D‘{z)  have  definitions  which  correspond 
exactly  to  Ii‘(z)  and  D‘(z). 

In  Eqs.  (2.34)  and  (2.35),  { cS,U).  c'ti(s)}  and  %>**,<?). 
(•’„(>),  are  si  ts  of  two  different,  but  consistent,  direct 
correlation  functions  with  corresponding  closures  (dis¬ 
cussed  m  the  next  section)  for  the  wall-particle  and 
particle-particle  interactions,  respectively.  The  func¬ 
tions  ,•,’,(,  )  and  <-{,(>•)  for  the  bulk  fluid  particles  are  ob¬ 
tained  in  an  independent  calculation  from 

(.’u<r)  -  <  i',(r) -pi/ij,  *  rJi  .  (2.30 

Iq i(r)  -  c) ](r)  -  K u p,  Iq j  *  I'u  .  (2.2< 

which  an  the  one-component  analogs  of  Eqs.  (2.28)  ..no 
(2.  29)  first  derived  in  a  seminal  paper  by  Werlheim. 
These  equations  carry  their  own  closures  for  //u(rl(r 
and  ru(rl(r  '■ Kn).  We  do  not  actually  solve  Eos.  (2.  36' 
and  (2.35)  since  nearly  exact  results  are  available  from 
the  work  of  Verlet  and  Weis6  for  hard  spheres  and  from 
the  study  of  Waisman.  Henderson,  and  Lebowitz1  for 
hard  spheres  against  a  wall.  Our  RLHNC  approximation 
(Sec.  Ill'  implies  that  /r2,(z)  is  the  exact  wall-particle 
total  correlation  function  for  hard  spheres  against  a 
hard  wall.  Equation  (2.37)  has  been  solved  by  Wei  - 
iheun.  m  the  mean  spherical  approximation. 6  w!ui< 
Paley  and  Ins  colleagues  have  treated  it  in  the  l.l!\C 


nr-  t  c  *  .r*(  *  lilt.  _>  IIIIVI  til.-1  o  iisnx  u  vatvvi  *t  m  nix 

h\i (>•)  -  e|i(r)  ■>  -  U~P-1  f  t/f  /  f  i/s  s  c'u(>l  .  and  QHNC  approximations. 6  using  the  Verlet-Wei 

'  tr-d  ......  ,U.  1...  ..a  ...  ll'r,  ,  O.r, 


(2.30) 

On  substituting  r  --  H:  ■  z .  /  -  H 2  -  v  and  on  taking  the  wall 
limit  R2—  x  .  we  find,  for  .-  •  0. 

l<hU)--ctilU)~2Kuzi>\  I  r/viV21(v)  /  r/ssr;,(<)  . 

-I J-vl 

(2.31) 

where  z  is  the  distance  of  the  center  of  the  dipolar  hard 
sphere  from  the  wall.  The  integral  between  the  limits 
-oc  and  00  may  be  simplified  by  noting  that  A2,(  v)  -  -  1 . 
when  -«<y<  0,  and  in  addition  w  -yl  -z  -  v.  when 
z  >0  and  -00  <  y  ■-  0.  On  changing  the  order  of  integra¬ 
tion  between  -*>  and  0. 

I  ilx  h\\(  v)  f  i/s  .<i'||(sl  I  <ls(z  -  s)  seals')  . 

a—  •'if-yl 

(2.32) 

Defining  the  functions 

fi*U)  =  f  c;,(s)  m/s  .  />*(<-)-  I  *c*„(s)s*rfs  .  (2.33) 

-'ll  -'ll 

the  wall-particle  ±  equations  become 

*«(■?■>  -t-St<^>  -  2irtfi,)>,i'e  !»*(*) -fl‘(c)  |-  [/'*(*■)  -/;*(••  !-:- 

•  2-K,,p;  f  ,/vfT21(\)|«*(^)-S*(;r-v  )|  . 

-'ll 

(2.341 

where  the  functions  B*(.  )  and  /)’(/)  are  entirely  deter¬ 


or\  for  the  hard-sphere  interactions.  We  are  then  ton 
lett  with  the  necessity  of  solving  only  Eq.  (2.34'  until  :• 
the  appropriate  closures  for  the  wall-particle  and  par¬ 
ticle-particle  direction  correlation  functions. 

The  constants  A’,,  and  A’2,  determine  the  dipole  mo¬ 
ment  w ,  and  the  electric  field  E2  through  the  relations 

trrlr-1  <<l.(2A'up1fl  Ji)  -  V-(~  ^ubt^ii'  <2. 


^ y-  '  A  -I  [  2<t',(2A'np, R  ,, )  4  t*f.(—  A,, p[  R  u )  ,  (2.39 

where  K  and  T  are  the  Boltzmann  constant  and  absolute 
temperature,  respectively,  and  <V..(Pi  R\{)  is  defined 
by: 

•y.-to,  A’  1 1 )  1  —  4 "Pi  I  rf,(r,  p,)  r*tlr  .  (2.401 

Equation  (2.38)  has  been  derived  by  Werthein  in  the 
mean  spherical  approximation,  but  its  extension  to  t:.' 
I.HN'C  tui  HEHNC)  approximation  is  straightforward: 
nevertheless  we  present  it  for  completeness  and  as  a 
prelude  to  the  derivation  of  (2.39).  which  Isbister  and 
1  reasier'  discussed  in  the  mean  spherical  approxm..: 
turn.  Tin'  derivation  of  Eq.  (2.38)  rests  on  the  .isvr:.;'- 
tod,  torn  of  e{>,(r)  (see  Eqs.  (3.1)  and  (3.6' 
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m2. 1)  .  .1.-  >•-  •  (2.  41) 

On  inserting  tins  m 

—3  /  .  (2.42) 

one  sees  that  rfjO')  is  short  ramanl  ami  tends  to  zero  as 
r  —  •*“ ; 

as  r  —  •*  .  (2.43) 

When  t ho  inverse  relation  n*f.  Kqs.  (2.  19)  and  ( 2 . 1 8)  | 


a  hie!:  ib  uirntie.d  U>  Kcj.  (2.39).  ll  should  In*  cmph.j- 
M/<  d  lh.it  I  'j  .  ( 2.  3b)  and  (2.39)  are  relations  drriv* 

I  nun  tin  a.-  >  inploi  ic  torins  of  rjj  and  in  tin*  Ms  an  : 
I.HNt  approximations.  The  dipole  moment  n:  ;  and  tn> 
ideeirif  held  are  determined  Iron)  An  and  A‘2I,  an* 
the  solutnms  t< j  hq .  (2.37)  have  been  obtained. 

The  constant  /i21  also  retrieves  Die  cuellicients  /  d. 
and  //21{.  ),  a  inch  appear  in  the  invariant  expansion  oj 
the  wall -particle  correlation  lunetions 

e;,  /»:,(. )  •  //?,(.-) m2,  n  -  .M2.  d  12.  *, 


<•?,(>•)  -  4  I  n,(-l>-</s  (2.44) 

)  a,, 

is  taken  to  the  limit  one  ..lsn  linds  that 


•Jh 

kT 


3  f 


<  ii(s)  s*tis 


--3A'U  j  [2(q,(s.  2/\,j,i.)  •  -A'uP,)  . 

(2.  45) 

where  we  have  used  the  analogs  of  t-qs.  (2.23)  and 
(2.24)  applied  to  e? ,(s)  in  the  last  stop.  The  result 
given  in  Eq.  (2.38)  follows  immediately  on  applying 
Eq.  (2.40). 


The  relation  between  A'.,  and  L ,  is  also  derived  from 
the  asymptotic  from  ot  (■?,(»•)  in  the  MS  and  LHNC  ap¬ 
proximations  [see  Eqs.  (3.1)  and  '3.8)  in  Sec.  Ill  j .  The 
equation  for  dipolar  mixtures  which  corresponds  to  Eq. 
(2.45)  is9 


Inn 

rjl-- 


»it»i  I  3A'2i  ( 

7x77  r  -  l‘n>  / 

r21'“  >  21  Vo 


[2<-;,(>)  -  rj,(s))ssrfs  . 


(2.46) 

where  we  have  not  canceled  the  r21  1!l  the  denominator, 
because  we  intend  to  take  the  wall  limit.  Substituting 
rzi  -  Rz  *  und  taking  the  limit  Rs  - Jf.  reduces  the  left 
hand  side  of  Eq.  (2.46)  to  ii/jE,,  K-T .  and  replaces  the 
upper  limit  of  the  integral  in  Eq.  (2.46)  by  * .  The 
Fourier  transform  of  Eq.  (2.29).  which  is  an  Ornstein- 
Zernike  equation  for  mixtures,  as  — 0,  yields 


?*„(*)  =  [i-a*„pI  ?;,(*)  |/1|,  a-) 


(2.47) 


from  which,  when  k-0,  we  have 


f  c‘z,(s)  s!ds  = 
•/o 


[i-a'„p:  f?,(o)i  | 


/i|,(s)s2rfs  . 


(2.48) 

In  the  above,  the  “'"represents  the  three  dimensional 
Fourier  transform 


/(/?)  =  fdrexpdk- r) /( i  r  j ) 


and  Mror  h. 

Since  /ij,(s)  is  a  short-ranged  function  equal  to  -  1  for 
s<R„,  it  readily  follows  that 

lim  4-  f  fr2l(s)  s2ds  -  -  3  .  (2.49) 

K  2"  *  y 21  -'o 

Hence. 


-  nji  (2  [  1  —  Aj|  P|  Ci|(0)  I  -  1 1  —  A' 11  Pi  Cj j(0)  |  | 


(2.50) 


by  inverting  Eq.  (2.23)  and  (2.24)  and  using  Eq.  (2.27.. 
Both  <Vj.(.  )  (for  the  closure  rules  considered  111  the  next 
section)  and  /<;,(.■)  are  short-ranged  functions  and  tenu 
to  zero  as  -  ■  .  and  the  asymptotic  form  of  //j,(.-)  1-, 
tlierefoi  c 

I1111  /it'd  I  3A,,  .  12.72 

winch,  follows  trom  Eq.  (2.27)  and  the  tact  that  l.';;1 
also  short  r.cmrd.  The  asymptotic  iorm  ol  the  wali- 
partiele  corn  latum  function  lor  the  It  LHNC  closure, 
discussed  in  tin  next  section,  coincides  with  that  ;  :  tn 
MS  approximation1  and  is  therolore  given  by 


Bin  .  E;.  <1,1  3A,,  D(2.  1) 

(2.  33' 

3o.  ,  h3lH2.  D 

v/  . 2^ 2A  1  i  / > j  1  j )  -  (Jl *“  A  u  /  Ajj)  1 

'2.  54' 

3FgD(2.  1)| 

■  r*  i , , p,  ft? *  (j| - X , .  j>, ft1,,!';  ’ 

(2.5'' 

where 

(2.  5t 

and 

h0R,A  lE;i  K\\ 

0  it,  j  (3  cosat>2  -  1)  '  ‘  tn  j 

(2. 57’ 

Since  the  electric  field  is  independent  of  the  distance 
from  the  wall.  htl{:.  E2.  ft,)  does  not  decrease  to  zero 
as  unless  the  magnitude  of  the  electric  field  or 

D(2,  1)  is  also  zero.  The  angular  average  ol  E:. 

is  however  zero  in  the  limit 

lim  /  h2 1(.'.  e2.  n,)r/n,  =  o 

e- -  J 

(2.58) 

since  the  angular  averages  of  £>(2,  1)  and  M2,  1) 
and 

are  zero 

lim  hzlU)  =  0  . 

(2.  59) 

We  shall  now  consider  the  details  of  these  equations  in 
the  context  of  the  mean  spherical  and  linearized  hyper - 
netted  chain  closure  rules. 

III.  THE  CLOSURE  RELATIONS 

The  closures  for  the  wall -particle  direct  correlation 
functions  arc  readily  derived  by  taking  the  wall  limit  ol 
the  corresponding  closures  in  the  bulk  fluid. 
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A.  MSA 

Tiu'  closure  tor  the  direr!  rrel.itum  function  is 


The  .drive  equations  (3.  121  were  solved  lUr.iliv-  i  . 
i  oiiij'  :ti  r  to  generate-  tlie  accompanying  .'lean  . 


<V(r)  ‘-’rir—f  />(1. 2)  .  tori  (3.1) 

K-l  y 

which  is  equivalent  to 

(•^(rl  0  .  for  i  •  /<„,  .  (3. 21 

This  is  unaffe  cted,  for  K:I.  on  taking  the  wall  limit 
of  C;,(  r) . 


B.  LHNC  approximation 

starting  with  the  liypernetted  chain  (UNO  approximu- 


rn„(2, 1) -lia42.  11  -  1111,^(2.  1*  -  fo6(2. 1) 

(3.3) 

and  using  the  invariant  expansions  of  ra6( 2,  1)  and 
h„A2.\)  .cf.  Eqs.  (2.131  ami  <2.14)1.  one  finds 

,„lr>  ^  -  .»k. 


a/' 


/ '( 2.  1)  . 


where  T‘s  is  the  spherically  symmetric  part  of  the  po¬ 
tential.  Expanding  the  logarithm  up  to  the  linear  term, 
and  collecting  and  comparing  coefficients  of  1.  D(2,  1), 
and  -M2,  1).  we  have 

r;„(r)  ■- /ii*(r)  -  tnjf^r)  -  V‘a(>)  kT  ,  (3.5) 


<-Sfi(r)^C(r)[l-^(r)“i‘7?^  ■ 


IV.  RESULTS  AND  DISCUSSION 

Tin  solution  of  Eq.  (2.51).  for  fi2,(.  .  E..  l.;'.  a.  - 
lam*  d  as  tile  confluence  of  three  distinct  compulation- 


(II  The  bulk  correlation  lunctions  wt  .  culcu:.a«-c. 
tin  maimer  described  by  Patey,  from  Eq.  (2.37  .  T  ■< 
reduced  dipole  moment  dependence  ul  this  relation  . 
'..lined  in  the  constant  A'**,  through  E<1.  (2.38k  i...r. 
ne  expressed  in  the  MSA  and  UL11NC  approximations  a- 


when-  r'n (/■<!*)  and  Af,(/fJ,)  are  the  values  of  I'i'dri  and 
h;.(r)  immediately  outside  contact.  The  value  o! 
was  adjusted  until  this  difference  assumed  the  desirid 
reduced  dijxde  moment.  The  iteration  ol  Eq.  (2.37 
was  performed  using  last  Fourier  transform  tecnniqu.  - 
and  mixed  solutions  to  speed  convergence.  The  bub. 

! iuici  correlation  functions  obtained  showed  txc< 
aureeme: i!  mfh  those  of  Patey." 


(2'  The  electric  field  which  emerges  Iron,  tta  w.d. 
was  determined  by  Eqs.  (2.9),  (2.30).  and  12.40'.  j  ,h. 
field  angle  a,  of  Fig.  1  was  taken,  in  separate  i-ulrul.. - 
lions,  as  O'.  45' ,  and  90".  The  field  angle  is  related 
to  the  wall-dipole  orientation  through  Eq.  (2. 101 .  the 
solution  of  which  appears  in  Fig.  1.  The  magnitude  ul 
the  reduced  electric  field 


E? ;  fin  E„(3 cos2fe? -  l)l/sR\ 


fj*(»-)=fcj,(r)[l  -Mldr)-':  • 

(3.7) 

Defining 

(3.8) 

b^<r)=//ZB(r)[l -frlilr)-'1'  . 

(3.9) 

with  b«B(r)  and  b'aB(r)  also  defined  by  equations  analogous 
to  Eqs.  (2.18),  (2.23),  and  (2.24),  respectively,  the 
last  two  equations  can  be  written  in  the  form 

clB(r)  =  b'lxB(r)  ,  r> RaB  .  (3.10) 

In  the  RLHNC,  Eq.  (3.5)  is  replaced  by  the  exact  closure 
for  hard  spheres  against  a  wall. 

In  the  waff  limit,  5?,(z)  =*£(z),  and  malting  use  of 
Eq.  (2.27),  the  closure  condition  for  the  wall -particle 
direct  correlation  function  becomes 

C|i(^)  =  [l  +*JiU)J[l  -^SiU)'*]  ,  *>0  •  (3.11) 

On  substituting  this  in  Eq.  (2.34),  we  have,  for  z>0, 

AfeUM* 1,W  -  1  h  tftiU)  2*K„P1  j*  [B‘(«)  -B‘U)  | 

'+[£>*<«) -£>*(z)  I  +  j  *|,U)[B*(«)-B*(|z-y|)Uv|  • 

(3. 12) 

When  g|,(<-)  - 1,  the  formal  equation  for  /i|[(z)  in  the 
MSA  approximation  is  recovered. 


was  taken,  as  in  the  earlier  work  of  Isbister  and 
F reusing  to  be  with  the  reduced  dipole  momen’ 
squared  mj8  fixed  at  either  0.5  or  1.0.  For  pur  post . 
of  comparison  of  these  parameters  to  a  molecular  sy-- 
tem,  the  reduced  dipole  moment  squared  for  HCi  (a  ; 

1.03  l»  at  275  "K.  assuming  a  diameter  of  3.5  A,  is 
approximately  0.65.  A  reduced  field  strength  o!  B  3 
for  this  system  corresponds  to  an  electric  field  of  near¬ 
ly  1.9  •  109  V'nt  or  a  surface  charge  density  a  of  1  elec¬ 
tronic  charge' 1000  A2. 

Wt  have  also  carried  out  calculations  at  the  same 
surface  charge  density  (or  electric  field  £2)  for  a  fluid 
at  a  reduced  density  of  0.  7  with  the  reduced  dipole  mo¬ 
ment  hi*  -  2.0.  These  numbers  correspond  approxi¬ 
mately  to  those  appropriate  for  liquid  water  (m,  -  1 .  B5 
D,  /?,,  =  2.  76  A)  at  room  temperature.  The  reduced 
electric  field  Bj  [which  contains  m,  and  in  its  defini¬ 
tion  (4.2)]  is  now  only  0.71.  (It  may  be  useful  for  the 
reader  to  bear  in  mind  that  it  is  an  artifact  of  our  defini¬ 
tion  of  E*  that  an  increase  in  the  dipole  moment  w,  re¬ 
sults  either  in  a  reduction  of  £j  when  the  surface  charge 
density  is  held  constant,  or  an  increase  in  the  surface 
charge  density  when  Ej  is  unchanged. ) 

(3)  The  spherically  symmetric  part  of  Eq.  (2.  51). 
t.  c. ,  lij,(.-),  was  determined  using  the  technique  of 
Waisman,  Henderson,  and  Lebowitz. 7  The  function 
/ij,(z)  was  first  obtained  from  Eq.  (2.28)  as  the  Percus- 
Yevick  solution  and  then  corrected  to  produce  an  essen- 
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ivy 


R.  5R, 


KKi.  The  wall- particle  distribution  functions  g  \*2>  E;.  1M 
as  a  function  of  the  distance  z  from  the  wall  for  different  orien¬ 
tations  of  the  fluid  dipoles  where  E*  =  s  3,  m*:  -  o. pi* 

-  U.373,  and  o.  -  h  .  k.  is  the  radius  of  a  fluid  dipole. - 

IU.HNV  approximat  ion, - MS  approximation. 


-2- 


7  ~ 


b  f\ 


5P, 


riU.  3.  Wall-particle  distribution  functions  for  the  system 
depicted  in  Fig.  if  except  fttst  a,  - 


FIG.  ).  Wall-particle  dislrihulion  functions  for  the  systems 
depicted  in  Figs.  2  and  ii  except  that  a1  =  ltti°. 


3R 


1IG.  f>.  'Hie  wall-particle  distribution  function  y  ■  .  E.,  li,i 

for  the  system  depicted  in  Fig.  I  except  that  ml  l.o. 


I  20  -  ’  ’  '  r~T 

I 


0  PC’  V- 


)  hi.  rhe  expansion  coefficients  ht  *z)  and  hl\  i<?>  as  a  June- 
lion  o|  the  wall-particle  distance  z  when  Ef  -  h  it  and  ,>T  -  o.  573, 
the  upper  and  lower  curves,  in  each  case,  are  lor  ».*  l.p 
and  0.5,  res|iectivelv. 
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l'iR, 


i'U,.  7 .  Hu*  w:  U- particle  i- 1 1  li-ic  inn  time  lion  e  .  i c,  E  ,  I.’,- 
as  a  hmclion  ol  the  distance  Ir--  ■  t-.c  wait.  tor  ditterent 
orientations  i.’,  o!  tin1  than  dipole-'  'a hen  £*  ‘‘.71.  m*  l.u, 

pf  0.7,  ami  0.  Che  rcd*iccd  parameters  -n*  ami  are 
those  appropriate  tor  liuuid  w.ner  ..[  room  tem|>erature  on, 

\ .  sr>  l),  Etl  .’.TO  \>.  the  reduced  electric  field  £*  0.71 

corresponds  to  a  surface  densilv  tit  1  eieetronie  charge  lnop  A  , 


tiully  exact  wall -particle  correlation  function  in  com¬ 
parison  to  Monte  Carlo  simulation.  The  reduced  density 
of  hard  spheres,  upon  which  hi,(r)  depends  solely,  was 
taken  as  p*  -piR,,  =0.  573. 

The  computations  were  combined  through  Eqs.  (2.23), 
(2.  24),  and  (3. 12)  to  yield  the  coefficients  of  the  angular 
functions  of  Eq.  (2.51).  The  integrals  over  bulk  corre¬ 
lation  functions  (2.33)  as  well  as  the  integral  of  Eq. 

(3.  12)  were  performed  by  trapezoidal  approximation. 

The  iterative  solution  to  Eq.  (3.12)  was  rapidly  conver¬ 
gent.  Finally,  the  expansion  coefficients  were  inserted 
into  Eq.  (2.51). 

In  Figs.  2-5  we  present  our  results  for  .g21(z.  E2,  ST ,) 
in  the  KLHNC  approximation.  The  expansion  coeffi¬ 
cients  hj,(r)  and  /i2l(r)  are  shown  In  Fig.  6.  A  compari¬ 
son  of  KLHNC  and  MSA  results,  for  a  field  angle  of  zero 
and  reduced  dipole  moment  squared  of  0.5  and  1.0,  are 
shown  in  Figs.  2  and  5,  respectively. 

In  the  case  of  m*2  =  0.  5.  F,'  -  8  3,  and  p*  -  0.  573 
tf21U,  E2,  Si,)  in  the  KLHNC  approximation  assumes 
negative  values  only  very  near  and  at  contact  for  a  di¬ 
pole  orientation  in  opposition  to  the  field.  For  higher 
values  of  w*2-  1.0,  at  the  same  reduced  electric  field 
and  fluid  density,  these  regions  are  more  pronounced 
and  extended  in  range  (Fig.  5).  An  increase  in  w*2 
from  0.5  to  1.0  at  constant  E*  and  p* ,  however,  corre¬ 
sponds  to  a  fourfold  increase  in  the  surface  charge  den¬ 
sity  from  1  electronic  charge  1000  A2  to  one  every 


250  A2.  Until  the  KLHNC  and  MSA  treatments  result 
m  distributions  showing  exhanccd  adsorption  lor  favor¬ 
able  di|xile  orientations,  but  file  repulsive  interaction.- 
of  tile  electric  field  with  unfavorably  aligned  dipoles  i- 
more  clearly  visible  in  the  KLHNC  approximation.  Tins 
theory,  like  the  MSA.  is  essentially  linear  in  character, 
and  cannot  prevent  the  distribution  functions  from  be¬ 
coming  negative  when  the  dipoles  are  aligned  against  li.< 
field,  it  at  the  same  time  the  theory  predicts  a  large  i  i  - 
hancemeni  ol  the  adsorption  of  dipoles  aligned  with  the 
field.  When,  for  example,  the  electric  field  is  perpen¬ 
dicular  lo  the  wall,  both  theories  predict  that  the  den¬ 
sity  proliles  (see  Figs.  2  and  5)  are  symmetrical  about 
the  protUe  tor  dipoles  perpendicular  to  the  held  f:  .. 

-  -  2  or  3;  21.  A  large  enhancement  of  dipoles  align- 1 
with  the  field  would  lead  to  an  equally  large  depletion  ,■ 
dipoles  oriented  against  the  field,  which  may  require 
negative  wall -particle  distribution  functions.  J-urthi  : 
improvements  beyond  the  RLHNC  approximation  mo.iI-; 
then  be  necessary.  Figure  7  shows,  however,  lii.it  tic 
KLHNC  theory  provides  plausible  density  prolhes  i  ven 
whim  the  square  of  the  dipole  moment  n.*‘  is  increased 
to  4.0  hut  the  surface  charge  density  is  maintainec.  at  1 
electronic  charge' 1000  A2. 
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